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ABSTRACT 

Context. The central kiloparsec of many local Luminous Infra-red Galaxies are known to host intense bursts of massive star formation, 
leading to numerous explosions of core-collapse supernovae (CCSNe). However, the dust-enshrouded regions where those supernovae 
explode hamper their detection at optical and near-infrared wavelengths. 

Aims. We investigate the nuclear region of the starbust galaxy IC 694 (=Arp 299-A) at radio wavelengths, aimed at discovering 



(— I ' recently exploded CCSNe, as well as to determine their rate of explosion, which carries crucial information on star formation rates, 

the initial mass function and starburst scenarios at work. 

I , Methods. We use the electronic European VLBI Network to image with milliarcsecond resolution the 5 .0 GHz compact radio emission 

O ' of the innermost nuclear region of IC 694. 

' Results. Our observations reveal the presence of a rich cluster of 26 compact radio emitting sources in the central 150 pc of the nuclear 

C/3 , starburst in IC 694. The large brightness temperatures observed for the compact sources indicate a non-thermal origin for the observed 

. radio emission, implying that most, if not all, of those sources are young radio supernovae (RSNe) and supernova remnants (SNRs). 

We find evidence for at least three relatively young, slowly-evolving, long-lasting RSNe (AO, A12, and A15) which appear to display 

unusual properties; suggesting that the conditions in the local circumstellar medium (CSM) play a significant role in determining the 

I radio behaviour of expanding SNe. Their radio luminosities are typical of normal RSNe, which result from the explosion of Type 

JiT . IlP/b and Type IlL SNe. All these results yield support for a recent (less than 10-15 Myr) instantaneous starburst in the innermost 

• regions of IC 694. 

in . 

0^ ' Key words. Galaxies: starbursts - luminosity function, mass function - individual: IC 694 - Stars: supernovae: general ~ Radiation 

' mechanisms: non-thermal - Radio continuum: stars 

o\ 
o 

^\ 1. Introduction Long Baseline Interferometry (VLBI) arrays, of the order of a 

. few milliarcsec at cm-wavelengths, is able to detect individual 

^ , The observed rate at which massive stars (M > 8 Mo) die as ^^^^^ supernovae at large distances in the local Universe. 
<^ ' CCSNe, vccsN, can be used as a direct measurement of the cur- 
rent star formation rate (SFR) in galaxies, and provides unique Starburst activity in the circumnuclear regions of (U)LIRGs 

information on the initial mass function (IMF) of massive stars, implies both the presence of a high number of massive stars and 

B While the rate at which stars die in normal galaxies is rather low ^ dense surrounding medium, so bright radio SNe are expected 

■ ■ ' (e.g., one SN is expected to explode in the Milky Way every ^ occur (|Chevalieij|1982t IChugaill 1997|), and high-resolution ra- 

-50 yr), the CCSN rate in Luminous and Ultraluminous Infra- ^io observations have shown that highly-extinguished CCSNe 

red Galaxies (LIRGs, L ,r = L[8 - 1000 Mm] > 10" L^; and do exist in the circumnuclear (r< 1 kpc) region of local 



ULIRGs, LiR > lO'^ L^: ISanders & Mirabell ([l996h ) is expected (U)LIRGs (Smith et al. 1998; Lonsdale et al. 2006; Colina et al 
to be at least one or tw o orders of magnitude larger than in nor- 200 1 ; Neffet al., 2QM; Perez-Torres et al, 2007; Kankare et al 



mal galaxies (ICondonlll992h . and hence detections of SNe in 2008|). Therefore, VLBI observations can set strong constraints 

(U)LIRGs offer a promising way of determining the current star on the properties of star formation m the dust-enshrouded envi- 

formation rate in nearby galaxies. ronments encountered in (U)LIRGs. 

However, the direct detection of CCSNe in the extreme den- Arp 299 consists of two interacting galaxies ( IC 694 and 

sities of the central few hundred pc of (U)LIRGs is extremely NGC 3690), which are in an early merger stage ("Keel & Wi 

difficult, since emission in the visual band suffers very signifi- 1995). At a luminosity distance of 44.8 Mpc (Fixsen et al. 199' 

cant extinction in those regions which contain large amounts of for Hq - 73 km s ' Mpc ', Arp 299 has a n infrared luminosity 

dust, and can at best yield only an upper limit to the true value of Lir ^ 6.7 x lO^L^r, (I Sanders et al"]|2003h . which is approach- 

I'ccsN- Fortunately, it is possible to directly probe the star form- ing the ULIRG category. The innermost ~150 pc nuclear region 

ing activity in the innermost regions of (U)LIRGs by means of of Arp 299-A (see Figure [1]) is heavily dust-enshrouded, thus 

high angular resolution (< 0.05 arcsec), high-sensitivity (<0.05 making the detections of SNe very challenging even at near- 

mJy) radio searches of CCSNe, as radio does not suffer from dust infrared wavelengths. Yet, Arp 299 hosts recent and intense star 

extinction, and the angular resolution yielded by current Very forming activity, as indicated by the relatively high frequency 
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of supernovae discovered at optical and near-infr ared wave- 
lengths in their outer, much less extinguished regions ("Forti et al.! 
1993; van Buren et al. 1994; Li et al. 1998; Yamaoka et al. 1998; 
biu et al.lll999l:lMattila et al.ll2005l) . 



The brightest component at infrared and radio wavelengths is 
IC 694 (A in the top panel of Figure[TJ hereafter Arp 299-A), ac- 
counting for ~50% of the total infrared luminosity of the system 
dAlonso-Herrero et al. 2000; Charmandaris et al. 2002), and for 
-70% of its 5 GHz radio emission dNeff et al. 2004). Numerous 
H II regions populate the system near star-forming regions, sug- 
gesting that star form ation has been occurring at a high rate 
for the last -10 Myr dAlonso-Herrero et al.l 12000). Given that 
IC 694 accounts for most of the infrared emission in Arp 299, 
it is t he region that is most likely to reveal new SNe dCondonl 
119921) . Since optical and near-infrared observations are likely to 
miss a significant fraction of CCSNe in the innermost regions 
of Arp 299-Al due t o large values of extinction [Ay ~ 34-40 
dGallais et al.ll2004l;lAlonso-Herrero et al.l2009i) l and the lack of 
the necessary angular resolution, radio observations of Arp 299- 
A at high angular resolution, high sensitivity are the only way 
of detecting new CCSNe and measuring directly and indepen- 
dently of models its CCSN and star formation rates. In fact. Very 
Long Baseline Array (VLBA) observations carried out during 
2002 and 2003 re sulted in the detection of five compact sources 
dNeff etal.ll2604h . one of which (AO) was identified as a young 
SN. 



2. eEVN observations and results 

We used the electronic European VLBI Network (e-EVN) 
dSzomoru..2006} to image Arp 299-A at a frequency of 5 GHz 
over 2 epochs, aimed at directly detecting recently exploded 
core-collapse supernovae through the variability of their com- 
pact radio emission (see appendix lAl for a detailed explanation 
of our observing strategy, calibration and imaging procedures, 
and source detection and techniques for flux density extraction). 
The attained off-source root-mean-square (r.m.s.) noise level was 
of 39 yuJy/beam and 24 yuJy/beam for the 8 April 2008 and 5 
December 2008 observations, respectively, and reveal the exis- 
tence of 26 compact components above 5 r.m.s. (see Figure [1]). 
Since the EVN radio image on 5 December 2008 is much deeper 
than the one obtained on 8 April 2008, it is not surprising that we 
detected a larger number of VLBI sources in our second epoch 
(25) than in our first one (15). This allowed us to go back to our 
first-epoch image and extract the flux density for the new compo- 
nents (A15 through A25 in Figure[T]i, and which show > 5rm.s. 
detections only in the December 2008 image). This procedure al- 
lowed us to recover four components above 3cr (A15, A18, A22, 
and A25), based on a positional coincidence with the peak of 
brightness of our second epoch of better than - 0.5 milliarcsec, 
i.e., much smaller than the synthesized interferometric beam). 

Our results demonstrate the existence of a very compact rich 
nuclear starburst in Arp 299-A and, in general , are in excellen t 
agreement with independent results reported bv lUlvestadl(l2009l) . 
The angular size encompassed by the radio emitting sources in 
Arp 299-A is less than 0.7"x0.4", corresponding to a projected 
linear size of (150x85) pc. To facilitate comparisons, we define 
here a fiducial supernova radio luminosity equal to three times 
the image r.m.s. in the 8 April 2008 epoch, which corresponds 
to 2.9x10^* erg s"' Hz"'. In this way, the radio luminosities for 
the VLBI components range between 1.1 (A25) and 7.3 (Al) 
and between 1.0 (A13) and 7.7 (Al) times the fiducial value, for 
the VLBI observations on 8 April 2008 and 5 December 2008, 
respectively (see Table[T]for details). 




Fig. 1. Top: 5 GHz VLA archival observations of Arp 299 on 24 
October 2000, displaying the five brightest knots of radio emis- 
sion in this merging galaxy. Middle and bottom: Contour maps 
drawn at five times the rm.s. of our 5 GHz eEVN observations 
of the central 500 light years of the Luminous Infrared galaxy 
Arp 299-A on 8 April 2008 and 5 December 2008, revealing a 
large population of relatively bright, compact, non-thermal emit- 
ting sources. The size of the FWHM synthesized interferometric 
beam was of (0.6 arcsec x 0.4 arcsec) for the VLA observations, 
and of (7.3 milliarcsec x 6.3 milliarcsec) and (8.6 milliarcsec x 
8.4 milliarcsec) for the EVN observations on 8 April 2008 and 5 
December 2008, respectively. To guide the reader's eye, we have 
shown in cyan the components detected only on the 5 December 
2008 epoch. 



3. Discussion 

The radio emission from the compact sources detected from our 
VLBI observations can be explained in principle within two dif- 
ferent physical scenarios: (i) thermal radio emission from super 
star clusters (SSCs) hosting large numbers of young, massive 
stars that ionize surrounding H II regions; (ii) non-thermal radio 
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Fig. 2. 5.0 GHz luminosity histogram for the VLBI components 
detected on 8 April 2008 (red) and 5 December 2008 (blue). The 
values centered at 0.0 correspond to the upper limits quoted in 
Table [T] See main text for details. 



emission from supernova remnants (SNRs) and/or young radio 
supernovae (RSNe), i.e., recently exploded core-collapse super- 
novae where the interaction of their ejecta with their surrounding 
circumstellar or interstellar medium (CSM or ISM, respectively) 
would give rise to significant amounts of synchrotron radio emis- 
sion. 

The existence of SSCs in Arp 299-A has been demon- 
strated by the ir apparent de tection using 2.2-yum adaptive op- 
tics imaging dLai et al] 1199 9). Further evidence in this direc- 
tion comes from Hubble Space Telescope (HST) FOC and 
NICMOS images, which reveal a popul ation of young stellar 
clusters in the central regions of Arp 299 ( Alonso-Herrero et alj 
[2OOO). The total 5.0 GHz radio luminosity in compact sources 
is 1.7x10^** erg s"' Hz"' and 2.0x10^*^ erg s"' Hz"' on 8 April 
2008 and 5 December 2008, respectively. However, the high 
flux densities in Tablefor most of the compact sources in Arp 
299-A ( <9 milliarcsec; see Table |2]i, indicate brightness tem- 
peratures, Tb, which largely exceed the thermal temperatures 
expected from SSCs (< 2 x 10'*K), thus ruling out a thermal 
origin for the compact radio emission traced by our eEVN ob- 
servations. 

Therefore, the observed radio emission must be due to young 
radio supernovae, SNRs, or both. We show in Table [1] the flux 
densities and luminosities for all the components identified in 
our observations, and classify the objects according to their vari- 
ability. The majority of sources do not show any evidence for 
significant variability which is consistent with their identifica- 
tion as SNR. Only three sources (A6, A 12, and A 15) show ap- 
preciable flux density variations between our two consecutive 
VLBI observations (see Table [T]for details), which is very dif- 
ficult to reconcile with their radio emission being due to SNRs. 
A6, A 12, and A15 are also detecte d in prev i ous V LBI obser- 
vations between 2003 and 2005 by lUlvestadI (1200^ . where he 
reports that he finds no evidence for significant variability at the 
~20% level. However, the individual flux densities at the various 
epochs of observation are not listed in his paper, which makes a 
detailed comparison with our data problematical. In fact, over 18 
months a recent RSN may have gone from very low to very high 
flux density values, and even to have faded away completely. 
Nevertheless, from his non-detection of A15 at 2.3 GHz, and our 
clear detection with an increasing flux density over two epochs 
at 5.0 GHz, we suggest that A15 is a relatively recent, slowly- 
evolving RSN. This behaviour is very similar to that displayed 



by AO (see below). Similarly, source A12 shows an increasing 
flux density at 5.0 GHz; and since it was also previously detected 
in 2003 and 2005, it too is likely to be a relatively r ecent and 
slowl y-evolving RSN. The nature of A6 is less clear. lUlvestadI 
(12009) detected A6 only at 2.3 GHz. This detection, together 
with the sudden drop in its 5.0 GHz flux density between April 
and December 2008 may be interpreted as the result of a CCSN 
leaving its young radio supernova phase (when its radio emis- 
sion is powered by interaction with the CSM) and now entering 
the supernova remnant phase (when it is powered by interaction 
with the ISM). We cannot exclude, however, the possibility that 
A6 is an X-ray binary or a microquasa r. In fact, A6 i s quite close 
to two X-ray sources reported by Zezas et al.) ((20031) (sources 14 
and 16 in their Table[T]). The combined absorption corrected X- 
ray luminosity of those sources is ~ 1 .9 x 10*' erg s"' . From the 
5.0 GHz radio luminosity of A6 in Table[Tl and assuming a spec- 
tral index of q'=-0.5 (typical of microquasars), the approximate 
total radio luminosity is ~ 4.8 x lO"'^ erg s"', resulting in a ratio 
of radio-to-X ray luminosity of ~ 2.5 x 10""^, which is a bit high, 
but still compatible with A6 being a microquasar. Approved new 
EVN high-sensitivity observations will allow us to confirm the 
nature of A6 in the near future. 

The maximum 5.0 GHz luminosities inferred by A12 and 
A15areof -1.6x10^^ erg s ' Hz ' and -7.6x10^^ erg s ' Hz ', 
respectively. These are typical values of radio emitting Type IIP, 
or Type lib SNe (Chevalier 2006). If confirmed, it would be the 
first time that such relatively faint radio SNe have been detected 
in the nuclear starburst of a local (U)LIRG. For comparison, Arp 
220 appears to contain essentially very b right radio supern ovae, 
which are identified with Type Iln SNe ( Parra et aP 2007). We 
note here that this is not simply a sensitivity issue, since some 
of the VLBI obser vations of Arp 220 had r.m.s. as low as 9//Jy/b 
(iParra et al.ll2007h . yet most of the objects detected there were 
identified, based on their large radio luminosities, as Type Iln 
SNe. We cannot exclude, however, that A12 and A15 have not 
yet reached their peak. In that case, their (peak) radio luminosity 
could be a few times higher, and given their slow evolution, these 
SNe could be Type IIL, or even Type Iln. 

In a previous paper, Neff et al. (2004) reported the detec- 
tion of five VLBI sources (AO to A4), within the central ~80 pc 
of Arp 299-A. We detect all of these components at 5.0 GHz. 
Components Al to A4 do not show significant variability in the 
eight months covered by our 5.0 GHz eEVN observations, which 
is consistent wit h them be i ng yo ung supernova remants, as is 
also suggested bv lUlvestadl(l2009h . 

AO was fi rst detected at 8.4 GHz, and identified with a radio 
supernova byl Neff et al.l (l2004t) . We have now detected it at 5.0 
GHz, more than five years after its discovery at 8.4 GHz. This 
implies that AO is a long-lasting, slowly-evolving, non-standard 
radio supernova, si nce most other examples evolve more rapidly 
dWeiler et alj|2002l) . Although rather uncommon, there are simi- 
lar cases reported i n the literature , both in normal galaxies, e.g., 
SN 1979C in MlOO lMontes et al I (I2000I). ari d in (U'lLIRGs. hke 
some of the RSNe in Arp 220 dParra et al1l2007l). In add ition, 
its non-detection at 2.3 GHz up to 2005 (lUlvestadI l2009h sug- 
gests the presence of a foreground absorber (e.g., a nearby H II 
region), as in the case of SN 2000ft in the LIRG NGC 7469 
(Alberdi et al. 2006; Perez-Torres et al. 2009). Finally, we note 
that A5 may also show similar properties, although its variabil- 
ity significance is lower. In short, the behaviour displayed by 
AO (and A5) whilst unusual, is not unknown for an RSN, and 
provides important information on how the interaction between 
the SN and CSM is proceeding, and thus probes the mass-loss 
history of the progenitor star. We therefore suggest that it is the 
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local CSM conditions which are primarily responsible for de- 
termining the power of the radio supernovae exploding in the 
nuclear starburst of Arp 299-A. 

e-EVN observations show that Arp 299-A hosts an ex- 
tremely prolific supernova factory, with radio luminosities typi- 
cal of Type lib. Hp, a nd IIL, and yielding s trong support to the 
scenario proposed by Alonso-Herrero et aTl (|2000), which sug- 
gests the existence of a recent (less than 10-15 Myr old) in- 
tense, instantaneous starburst. We find evidence for at least three 
slowly-evolving, long-lasting, non-standarnd RSNe (AO, A 12, 
and A15), which is very suggestive of the local CSM condi- 
tions playing a main role in determining the radio behaviour 
of the exploding SNe. Our current monitoring of Arp 299-A 
with the eEVN at 5.0 GHz, which is scheduled to continue until 
the end of 2010- should allow us to detect any new radio su- 
pernova, and therefore test whether the IMF in Arp 2 99-A is 
top-h eavy, in contr ast with n ormally assumed Salpeter (Salpeter 
Il955h . or Kroupa jKroupal 12001.) IMFs, where the production 
of massive stars (M > 8 M©) that eventually result in CCSNe 
is low compared to the production of less massive stars. There 
seems to be evidence tha t this might also be the ca se of M 82 
dPoane & Mathewslll993l) and Arp 220 ( Parraet al .'2007), and 
theoretically it is expected that in the warm dense ISM condi- 
tions within a (U)LIRG, t he IMF should inde ed be top-heavy 
due to a larger Jeans mass dKIessen et al.ll2007[) . 
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Table 1. Compact radio emitting sources in the central region of Arp 299-A 



Source 


Source 


Act" 


A6 


s„ 




U/IO^" erj 


; S-' Hz-' 




Name'' 


Type 


(J2000.0) 


(J2000.0) 


8 Apr 2008 


5 Dec 2008 


8 Apr 2008 


5 Dec 2008 




AO 


SN 


33.6212 


46.707 


318±42 


446+33 


7.9+1.0 


11.1+0.8 


2.4 


Al 


SNR 


33.6199 


46.699 


855±58 


901+51 


21.3+1.4 


22.4+1.3 


0.6 


A2 


SNR 


33.6219 


46.655 


708±53 


713+43 


17.6+1.3 


17.7+1.1 


0.1 


A3 


SNR 


33.5942 


46.560 


398±44 


353+30 


9.9+1.1 


8.8+0.7 


0.8 


A4 


SNR 


33.6501 


46.537 


558±48 


628+40 


13.9+1.2 


15.6+1.0 


1.1 


A5 


SN? 


33.6176 


46.693 


278±41 


143+25 


6.9+1.0 


3.6+0.6 


2.8 


A6 


SN 


33.6206 


46.597 


208±40 


<72 


5.2+1.0 


<1.8 


3.4 


A7 


SNR 


33.6300 


46.786 


496±46 


468+34 


12.3+1.2 


11.6+0.8 


0.5 


A8 


SNR 


33.6306 


46.401 


226±41 


264+27 


5.6+1.0 


6.6+0.7 


0.8 


A9 


SNR 


33.6306 


46.620 


294±42 


282+28 


7.3+1.0 


7.0+0.7 


0.2 


AlO 


SNR 


33.6392 


46.551 


550±48 


436+32 


13.7+1.2 


10.9+0.8 


2.0 


All 


SNR 


33.6403 


46.583 


300±42 


351+30 


7.5+1.0 


8.7+0.7 


1.0 


A12 


SN 


33.6495 


46.590 


449±45 


639+40 


11.2+1.1 


15.9+1.0 


3.2 


A13 


SN? 


33.6531 


46.733 


251±41 


118+25 


6.2+1.0 


2.9+0.6 


2.8 


A14 


SNR 


33.6825 


46.571 


292±42 


260+27 


7.3+1.0 


6.5+0.7 


0.6 


A15 


SN 


33.5991 


46.638 


159±40 


304+28 


4.0+1.0 


7.6+0.7 


3.0 


A16 


unci. 


33.6149 


46.516 


<117 


147+25 


<2.9 


3.7+0.6 


1.2 


A17 


unci. 


33.6179 


46.561 


<117 


179+26 


<2.9 


4.5+0.6 


2.4 


A18 


unci. 


33.6192 


46.464 


151+40 


129+25 


3.8+1.0 


3.2+0.6 


0.5 


A19 


SN? 


33.6196 


46.659 


<117 


191+26 


<2.9 


4.8+0.6 


2.8 


A20 


unci. 


33.6278 


46.789 


<117 


146+25 


<2.9 


3.6+0.6 


1.2 


A21 


unci. 


33.6291 


46.836 


<117 


133+25 


<2.9 


3.3+0.6 


0.6 


A22 


unci. 


33.6354 


46.678 


173+40 


217+26 


4.3+1.0 


5.4+0.7 


0.9 


A23 


unci. 


33.6360 


46.560 


<117 


137+25 


<2.9 


3.4+0.6 


0.8 


A24 


unci. 


33.6361 


46.764 


<117 


166+25 


<2.9 


4.1+0.6 


2.0 


A25 


SN? 


33.6524 


46.701 


132+40 


209+26 


3.3+1.0 


5.2+0.7 


1.6 



" Coordinates are given with respect to a(J2000.0)= 11:28:00.0000 and S (J2000.0) = 58:33:00.000, and were obtained from the 5 December 
2008 image. The positions for those sources also detected in 8 April 2008 coincide within the errors ( < 0.5 mas) with all of them. 
The uncertainty in the reported flux density for the detected compact components corresponds to Icr, where cr was determined by adding in 
quadrature the off-source r.m.s in each image and a 5% of the local maximum, to conservatively account for possible inaccuracies of the eEVN 
calibration. 

We define the significance of the flux density variability between the two consecutive epochs as V =| 5 Dec - Apr I / -^''"dcc ''"Ipr' where 5 Apr 
and cTApr (5 Dec and trDec) are the values in columns 5 and 6 (7 and 8), respectively. 

Source names are given in right ascension order, except for the five components reported previously (AO through A4) by jNeff et al.ll2004h . 



M.A. Perez-Torres et al.: An extremely prolific supernova factory in the buried nucleus of the starburst galaxy IC 694, Online Material p 1 



Online Material 



M.A. Perez-Torres et al.: An extremely prolific supernova factory in the buried nucleus of the starburst galaxy IC 694, Online Material p 2 



Appendix A: eEVN observations of Arp 299-A 

A.1. Observing strategy 

We observed the central regions of Arp 299-A at a frequency of 
5 GHz in two epochs using the EVN. Our first epoch on 8-9 April 
2008 (2008.99; experiment code RP009) included the follow- 
ing six antennas (acronym, diameter, location): Cambridge (CM, 
32 m. United Kingdom), Medicina (MC, 32 m, Italy), Jodrell 
Bank (JB, 76 m. United Kingdom), Onsala (ON, 25 m, Sweden), 
Torun (TO, 32 m, Poland), and Westerbork array (WB, 25 m. 
The Netherlands). Our second observing epoch on 5 December 
2008 (2008.340; RP014A experiment) included, in addition the 
EVN antennas Effelsberg (EE, 100 m, Germany), Knockin (KN, 
25 m. United Kingdom), and Shanghai (SH, 25 m, China). 

Both observing epochs consisted of e-VLBI phase- 
referenced experiments, using a data recording rate of 5 12 Mbps 
with two-bit sampling, for a total bandwidth of 64 MHz. The 
data were correlated at the EVN MklV Data Processor at JIVE 
using an averaging time of 1 s. First epoch observations con- 
sisted of ~8.0 hr on target. The telescope systems recorded both 
right-hand and left-hand circular polarization (LCP and RCP) 
which, after correlation, were combined to obtain the total inten- 
sity mages presented in this paper. 4.5 minute scans of our target 
source, Arp 299-A, were alternated with 2 minute scans of our 
phase reference source, Jl 128+5925. 3C345 and 4C39.25 were 
used as fringe finders and band-pass calibrators. Our second 
epoch consisted of ~4.5 hr on target. The telescope systems also 
recorded in dual polarization, and 4.5 minute on-source scans 
were alternated with 1 minute scans of Jl 128+5925. The bright 
sources 3C84, 3C138, 4C39.25 and 3C286 were used as fringe 
finders and band-pass calibrators. (We note that the inclusion of 
EE in the second observing epoch allowed us to get a much bet- 
ter r.m.s., and this in spite of the significantly less amount of total 
observing time.) 

A.2. Data calibration and imaging 

We analyzed the correlated data for each epoch using 
the NRAO Astronomical Image Processing System {AIPS; 
http://www.aips.nrao.edu). The visibility amplitudes were cali- 
brated using the system temperature and gain information pro- 
vided for each telescope. Standard inspection and editing of the 
data were done within AIPS. The bandpasses were corrected us- 
ing the bright calibrator 4C39.25. We applied standard correc- 
tions to the phases of the sources in our experiment, including 
ionosphere corrections (using total electron content measure- 
ments publicly available). 

Due to the limited bandwidth of CM and KN, the usable 
data of these antennas was found in a single sub-band with very 
noisy edges. To improve the quality of the bandpass calibra- 
tion, we removed the edges of CM for the first epoch, and the 
edges of CM and KN for the second epoch. The instrumental 
phase and delay offsets among the 8-MHz baseband converters 
in each antenna were corrected using a phase calibration deter- 
mined from observations of 4C 39.25. The data for the calibra- 
tor J 1128+5925 were then fringe-fitted in a standard manner 
We then exported the J l 128+5925 data into the Caltech imag- 
ing program DIEMAP (IShepherdetalJll995h for mapping pur- 
poses. We thus determined gain correction factors for each an- 
tenna. J 1 128+5925 showed a flux of ~ 0.38 Jy at 5.0 GHz at both 
epochs, and displayed a point-like structure at mas scales. After 
this procedure was completed within DIEMAP, the data were 



read back into AIPS, where the gain corrections determined by 
DIEMAP were applied to the data. 

The final source model obtained for J 1128+5925 was then 
included as an input model in a new fringe-fitting search for 
Jl 128+5925, thus removing the structural phase contribution to 
the the solutions of the delay and fringe rate for our target source, 
Arp 299-A, prior to obtaining the final eEVN images shown in 
Eigure [1] The phases, delays, and delay-rates determined for 
J 1128+5925 were then interpolated and applied to the source 
Arp 299-A. This procedure allowed us to obtain the maximum 
possible accuracy in the positions reported for the compact com- 
ponents in Eigure [T] We note, though, that this has no impact 
in the final images, since Jl 128+5925 is essentially point-like at 
the angular resolution ( < 9 milliarcseconds) provided by our 5.0 
GHz e-EVN observations, and therefore its phase-contribution is 
negligible. 

The imaging and deconvolution of Arp 299-A was done with 
the AIPS task IMAGR (see EiguresfHand lA.ll i. using a natural 
weighting scheme and a ROBUST parameter equal to zero. This 
scheme was a good compromise to maximize a very high an- 
gular resolution and a very low rm.s. values in our images. We 
point out that we kept the averaging integration time to 1 s and 
used a maximum channel bandwidth in the imaging process of 
8 MHz, which results in a maximum degradation of the peak re- 
sponse for the component furthest away from the phase center of 
less than 5% a nd prevents artificial smearing of the images (e.g., 
iBridle & Schw ab 1999). 

No self-calibration of the phases, nor of the amplitudes, 
was performed on the data, since the peaks of emission were 
very faint for such procedures. We emphasize here two pow- 
erful aspects of the use of the phase-reference technique: (i) it 
effectively increases the integration time on a source from min- 
utes to hours, thus increasing significantly the array sensitivity, 
and allowing the detectio n and imaging of very faint objects 
(iBeaslev & Conwavll 19951) and (ii) it retains the positional infor- 
mation of the compact components in Arp 299-A with respect to 
the phase-reference source Jl 128+5925, which lies 0.86° from 
the target source. 

A.3. Source detection, identification, and flux density 
extraction 

Our 5.0 GHz eEVN observations on 8 April 2008 and 5 
December 2008 resulted in off-source r.m.s noises of 39//Jy/b 
and 24/iJy/b, respectively. The difference in sensitivity was pri- 
marily due to the addition of the Effelsberg antenna in our sec- 
ond observing epoch, but also to the fact that there were more 
antennas than in our first epoch. We considered as real sources 
those individual peaks which had a signal-to-noise ratio equal 
to, or above five times the r.m.s. noise in either of our two ob- 
serving epochs. Since our second observing epoch was much 
more sensitive than the first one, this resulted in a larger num- 
ber of detections of compact components (26) with respect to 
the first epoch (15). We then put small boxes in the 8 April 2008 
image, around the positions of those components detected on 5 
December 2008, to extract the flux density of the counterparts 
initially detected only in the second epoch (components A15 
through A25 in Eigure[TJ. This procedure allowed us to recover 
four components above 3cr (A15, A18, A22, and A25), based on 
a positional coincidence with the peak of emission in our second 
epoch to better than 0.5 milliarcsec; and for those components 
not detected on 8 April, to derive 3cr upper limits to the emis- 
sion. 
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Fig.A.l. 5.0 GHz eEVN observations of the central ~150 pc 
years of the Luminous Infrared galaxy Arp 299-A on 8 April 
2008 ( top ) and 5 December 2008 ( bottom ), which reveals a large 
population of relatively bright, compact, non-thermal emitting 
sources. The colour scale goes from -SOjxJy/b up to 250 fiJy/b. 
The size of the FWHM synthesized interferometric beam was of 
of 7.3 mas X 6.3 mas and 8.6 mas X 8.4 mas, respectivetly. 

The flux density extraction was carried out within AIPS, us- 
ing task IMFIT, which fits Gaussian components to the image. 
Namely, we put small boxes (the few inner pixels of each source) 
around each local peak of brightness (above five times the rm.s. 
noise) and fitted single Gaussian components to each of them. 
Since the fitted integrated flux densities did not differ signifi- 
cantly from the peaks of emission (less than 1 %), it is clear that 
there is no evidence that any of the compact components are re- 
solved even with the very high angular resolution of a few mil- 
liarcsec that the eEVN delivers. The total flux densities which 
are listed in Table[T]are thus the peak flux densities found in our 
images. 
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Table 2. Brightness temperatures'' of the VLBI sources in Arp 299-A 



Source 




8 April 2008 






5 December 2008 




Name 


S. 0"Jy) 


a(mas) 


fo(mas) 




Sr (pJy) 


a(mas) 


fo(mas) 




AO 


318±42 


< 2.8 






446+33 


5.3 


3.0 


1.8x10* 


Al 


855±58 


< 3.6 






901+51 


2.1 


1.8 


1.6x10' 


A2 


708±53 


< 2.5 


< 1.5 


> 1.3 X 10' 


713+43 


2.1 


1.6 


1.4x10' 


A3 


398±44 


< 2.2 






353+30 


< 6.5 


< 0.3 


> 1.2 x 10' 


A4 


558±48 


< 3.6 


< 2.3 


> 1.0 X 10 


628+40 


< 2.8 


<0.5 


> 3.0 X 10' 


A5 


278±41 


< 11.1 


< 1.2 


> 1.3 X 10* 


143+25 








A6 


208+40 


< 7.6 


< 1.4 


> 1.3 X 10* 


<72 






••• 


A7 


496+46 


< 2.5 


< 3.3 


> 4.0 X 10* 


468+34 


3.4 


2.0 


4.6 X 10* 


A8 


226+41 


< 4.8 






264+27 


< 3.3 


< 0.7 


> 7.6 X 10' 


A9 


294±42 


< 3.6 


<4.4 


> 1.2 X 10* 


282+28 








AlO 


550±48 


< 3.2 


< 3.9 


> 2.9 X 10* 


436+32 








All 


300+42 


< 5.5 






351+30 


<4.5 


< 1.8 


> 2.0 X 10 


A12 


449+45 


< 3.1 


< 3.7 


> 2.6 X 10 


639+40 


2.3 


1.5 


1.2x10' 


A13 


251+41 


< 3.3 


<4.2 


> 1.2 X 10* 


118+25 


< 5.3 






A14 


292+42 


< 5.0 


< 3.7 


> 1.1 X 10 


260+27 


< 4.9 






A15 


159+40 








304+28 


< 2.8 


< 1.6 


> 4.5 X 10° 


A16 


<117 








147+25 


5.6 


3.2 


5.5 X 10' 


A17 


<117 








179+26 


7.3 


3.9 


4.2 X 10' 




1 ^1 -i-AC\ 








iZV + Zj 


< 5.1 


<5.5 


> D.l X lU 


A19 


<117 








191+26 


< 4.0 


< 2.0 


> 1.6 X 10* 


A20 


<117 








146+25 


< 7.9 


< 3.5 


> 3.5 X 10' 


A21 


<117 








133+25 


< 7.4 


< 2.9 


> 4.1 X 10' 


A22 


173+40 








217+26 


< 8.6 


<4.3 


> 3.9 X 10' 


A23 


<117 








137+25 


< 8.8 


< 2.9 


> 3.6 X 10' 


A24 


<117 








166+25 


< 8.6 


< 7.0 


> 1.8 X 10' 


A25 


132+40 








209+26 


< 8.6 


<4.4 


> 3.7 X 10' 



The brightness temperatures shown for the 5.0 GHz VLBI source components in Arp 299-A were calculated using the flux densities in 
Table [T]and the angular sizes quoted here. We derived the brightness temperatures from the general formula: = (2c^ /k) ByV^^ , where 
By is the intensity, in erg s"' Hz"' str"'. Since By depends on the measured flux density, Sy, and on the deconvolved angular size of each 
VLBI component (obtained by fitting them to elliptical Gaussians, characterized by their major and minor semi-axis, a and b). Therefore, the 
above formula can be rewritten as Jb = (2 c'/k) By v^- = 1.66 x 10' SyV^-(a by\ with 5v in mjy, v in GHz, and a and b in milliarcseconds, 
respectively. 



